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We derive a closed-form solution for the Green's function for the wave equation of a static (with 
respect to an undragged, static observer at infinity) scalar charge in the Kerr space-time. We 
employ our solution to obtain an analytic expression for the self-force on such a charge, comparing 
our results to those of Ref. []]. 



I. INTRODUCTION 

Understanding the Green's function on the Kerr ge- 
ometry is important in a number of fundamental the- 
oretical problems such as the radiation reaction prob- 
lem [5] 0], miniature black hole creation [3] and in 
the Kerr/CFT correspondence 5~:7 a - The wave equa- 
tions satisfied by Green's functions for black hole space- 
times are extremely complicated, even in the most sim- 
ple case of a Schwarzschild black hole. However, in 
the Schwarzschild case, despite not having closed-form 
solutions for the Green's functions, the static or zero- 
frequency scalar Green's function is known in closed form 
and may be used to obtain, for example, analytic expres- 
sions for the vacuum polarization on the horizon of the 
black hole [5] or the self-force on a static scalar charge 

For the Kerr space-time, as in the Schwarzschild 
case, the wave equation is fully separable [10] and we 
may therefore obtain a mode-sum representation for the 
Green's function in terms of spheroidal functions and 
complicated radial functions that must be computed nu- 
merically. The zero-frequency mode simplifies signifi- 
cantly though still a highly non-trivial double mode-sum 
over a product of four associated Legendre functions, not 
all of integer order. The only closed-form expression for 
a Green's function was obtained by Linet |llj . where he 
derived the scalar Green's function with one point on the 
pole of the Kerr black hole and Frolov [T2] used this re- 
sult to calculate analytically the vacuum polarization on 
the pole of the Kerr black hole. 

In this paper, we shall obtain the closed-form represen- 
tation of the Green's function for a static scalar charge in 
the Kerr space-time, i.e., a scalar charge at fixed spatial 
Boyer-Lindquist coordinates. This derivation relies on 
two interesting results concerning the associated Legen- 
dre functions, one of which we prove in this paper and the 
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other we prove in Ref. [12] ■ We validate our expression 
for the Green's function by calculating the self-force on 
a static scalar charge in the Kerr space-time, obtaining 
the result conjectured by Burko and Liu [T]. Therefore, 
this paper is also the first proof of the result claimed in 
Ref. [I]. 



II. THE KERR GREEN'S FUNCTION 



The Kerr metric in Boyer-Lindquist coordinates is 
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The event and Cauchy horizons arc defined by vanishing 
A, which has two roots r± — M ± \J M 2 — a 2 , the pos- 
itive subscript being the event horizon and the negative 
being the Cauchy horizon. Moreover, the Killing vec- 
tor d/dt is the null generator of the static limit surface, 
given by r e = AI + \JM 2 — a 2 cos 2 8. Inside this surface, 
the ergosphere region, d/dt becomes space-like and so 
local radial null cones tip over forcing particle's to move 
in the direction of the black hole's rotation, i.e., a par- 
ticle inside the static limit surface cannot remain static. 
We would therefore expect our closed-form expression for 
the static scalar Green's function to diverge as the static 
limit surface is approached. 

The Green's function for a scalar particle in the Kerr 
space-time satisfies the following wave equation: 
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The wave equation on this metric admits a fully separa- 
ble solution [10] and so the Green's function may be writ- 
ten as the following mode-sum representation in terms of 
separated mode-solutions: 
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where S UJ i m (cos0) are the normalized spheroidal func- 
tions satisfying 
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Eq. (2.5 1 and the normalization condition (2.6) imply 



that for u> = 0, the spheroidal functions rec 



uce to the 



normalized associated Legendre polynomials, i.e., 
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The radial part of the Green's function, Xcoim, satisfies 
the inhomogeneous equation 
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For a static source, the derivatives with respect to t van- 
ish in Eq. (2.3). Hence the static Green's function satis- 
fies 
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In this case, the mode-sum expression reduces to the 
zero- frequency mode of Eq. ( 2.4 1 (modulo a factor of 2ir) 



resulting in a simplified mode-sum solution in terms of 
associated Legendre functions 
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where Xim( r , r ') satisfies 



{dr ( 



d / . d 
A dr 



A 



l(l + l)}xim(r,r') = -6(r-r'). 

(2.11) 



Changing the radial variable 
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where b = \J M 2 — a 2 , 



we may rewrite the radial equation as 
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where 7 = a/b. We note that the regular singular point 
has now been shifted from r = r+ to 77 = 1. The solu- 
tions of the corresponding homogeneous equation are the 
associated Legendre functions of pure imaginary, non- 
integer order, -P; ±J7m (?7) and Qf nm (r]). The boundary 
conditions imposed correspond to different Green's func- 
tions which in turn correspond to different combinations 
of these associated Legendre functions. For the retarded 
Green's function, for example, one typically chooses in- 
going radiation boundary conditions at the horizon and 
outgoing radiation boundary conditions at infinity. How- 
ever, since we are considering a static particle, we cannot 
impose radiation boundary conditions. Rather, for the 
retarded Green's function, we require that the inner so- 
lution be regular on the future event horizon in ingoing 
Kerr-Newman coordinates while the advanced Green's 
function corresponds to the inner solution being regular 
on the past event horizon in outgoing Kerr-Newman co- 
ordinates. Near the horizon, the two candidates for the 
inner solution, P^ nm (rj), are oscillatory 
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In ingoing Kerr-Newman coordinates (r, 0, </>), where <f> 
is related to Boyer-Lindquist (f> by 
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as 77 — > 1 + . It is clear from these asymptotic expressions, 
that the first solution oscillates infinitely fast in ingoing 
Kerr-Newman coordinates whereas the second solution is 
regular. The appropriate inner solution for the retarded 
Green's function is therefore P; _ "™ 7 (?7). Similarly, for 
the advanced Green's function, the appropriate choice is 

(7). 

Having fixed our choice for the inner solution, the 
choice for the outer solution is determined by regular- 
ity at infinity; we may choose either Q\ mi (rf) or QJ lmi 
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since both are regular at infinity and differ only by a 
constant multiplier for each I and m. This constant mul- 
tiplier will be accounted for by the normalization. The 
solution of the retarded inhomogeneous radial equation 
Eq.(2.13) may be written formally as 
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where N = -(z 2 - l)W[Pf(z), Q((z)] and N = -(z 2 - 
l)W[Pf v (z) : Qj v (z)\ are constant and W[P,Q] is the 
Wronskian. For the pair of solutions {Pr"(v);QUv)}, 
the Wrosnkian takes a particularly simple form [14] : 
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The final form of the static retarded Green's function is 
therefore 
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where Acf> = </>' — <f>. 

The form of this Green's function looks similar to those 
that are known in closed-form. In particular, consider- 
ing the three-dimensional Euclidean Green's function in 
spheroidal coordinates gives us a product of four Legen- 
dre functions of integer order similar to the mode-sum 
above. In |15j . we have shown how one may also ob- 
tain a quasi-closed form for a product of four Legendre 
functions of non-integer order and degree by considering 
the Green's function on a three-dimensional cosmic string 
space-time in spheroidal coordinates. In the static Kerr 
case however, one has both integer and non-integer (pure 
imaginary) order Legendre functions and it is precisely 
this mixing that has made this problem so elusive for so 
long. In the following sections we shall show that the 
mode-sum above may indeed be performed leading to 
a completely closed-form expression for the static Kerr 
Green's function. In this paper, we offer one particu- 
larly neat application of this closed-form Green's func- 
tion, where we obtain an analytic expression for the self- 
force on a static scalar charge in the Kerr space-time and 
we compare our result to the conjectured result of Ref. [T] . 
We note also that the result of Ref. [1] relies on a con- 
jecture for the MSRP parameter which has recently 
been computed, along with higher order parameters, for 
the Schwarzschild [16] and Kerr [17] space-times. 
III. THE STATIC KERR GREEN'S FUNCTION 
IN CLOSED FORM 

In this section, we derive the static Kerr Green's func- 
tion in closed form for the field satisfying retarded bound- 



ary conditions; the derivation for the advanced field is 
analogous. Obtaining the closed-form solution relies cru- 
cially on two formulae involving the associated Legendre 
functions, results which may well prove useful in other 
contexts. 

In Ref. [15] , it is shown that a product of Associated 
Legendre functions of arbitrary order may be written in 
the following integral form: 
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This is a generalization of a result obtained in Ref. [TS] 
for a product of Associated Legendre functions of in- 
teger order. The derivation in [13j involves obtaining 
a closed-form solution for the wave-equation on a di- 
mensionally reduced Bertotti-Robinson space-time and 
equating with the equivalent mode-sum expression. Em- 



function (2.201, we obtain 



ploying our result (3.1 1 in our expression for the Green's 
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where we have changed the order of summation in the 
expression above. 

In the Appendix, we derive the following summation 
formula for the product of three Legendre functions 
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and x = cos A must lie in the range 6 — 6' < \ < min{6> + 
9',2tt — 9 — 9'}. For A outside of this range, the series 
sums to zero. An equivalent form of this result can be 
found in Ref. [19] though our derivation generalizes the 
method of Baranov [20 j who proves the m = case. This 
summation formula allows us to perform the Z-sum in our 
Green's function expression yielding 
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where the integral vanishes whenever x lies outside the range (x-,x+) where x± — cos 6 cos 0' ± sin sin 0'. The 
appearance of the integral may be greatly improved with the introduction of the independent variable defined by 
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and we have written the cosine in the numerator in Eq.(3.6) as a sum of exponentials. It is easy to see that these 
integrals are invariant under the change ^ — > — \& in the integrand so we obtain the single integral 

G static (x,x0 = 4, I" t exp ( im ^ Z ^ - 7 cosh' 1 ( X ))) 
The m = term gives the static, axisymmetric Green's function 
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which has been previously noted by Linet |21j using axisymmetric potential theory. The m-sum can be performed 
using the Fourier representation of a periodic delta function 
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The integral is now somewhat trivial, since we pick up contributions only at the zeros of the delta function, though this 
is complicated by the fact that the zeros are solutions of a transcendental equation, which we cannot solve analytically. 
Nevertheless, the integral may be performed using the definition of the delta composition: 
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where the sum is over all roots xi of f(x) = 0. In particular, we have 
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where cos Xi = cos 9 cos 9' + sin 9 sin 8' cos and are solutions of the transcendental equation 
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It turns out to be useful to re- write Eq.(3.16) in the form 
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(t? 2 + r/ 2 - 2 W ' cos A, - sin 2 A,) 1/2 = W - l) 1/2 (r/ 2 - l) 1 / 2 sinh - A0|) (3.19) 

where we have assumed, without loss of generality, that a > 0. This expression may now be used to greatly simplify 
our Green's function: 
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so that the are roots of /(^). Solving for the turning 
points, f'(^f) = is equivalent to solving 
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For points outside the ergosphere, which are the only 
physically relevant points to consider for a static Green's 
function, we have the following inequalities 

A 1/2 > a sin 9, A' 1/2 > a sin 9' , for a > 0, 

<9 < 7T, (3.23) 



so that the turning points $ satisfies the inequality 
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Since sinh a; is monotonic, there is precisely one turning 
point in the interval 
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Furthermore, it is easy to verify that this turning point 
is a minimum. Hence there are, in general, two solu- 
tions to the transcendental equation (3.16) which are 



approximately symmetric about the turning point. It 
turns out, however, that one of these solutions ought to 
be discarded since its inclusion in the Green's function 
yields the wrong Hadamard singularity structure. We 
take the solution that lies to the right of the turning 
point (which means we can discard the absolute value 
sign enclosing the denominator of Eq.(3.20)) since only 
this root yields the correct Hadamard singularity struc- 
ture for the retarded field. We can justify this mathe- 



matically by virtue of the ambiguity in the integral ( 3.10 ) 



which may be taken over any interval of length 2ir since 
the integrand is periodic; thus we simply choose any in- 
terval which includes the root of Eq. (3.161 which gives the 



correct singularities but which excludes the other root. 
Taken over such an interval, the integral in Eq. (3.131 only 



picks up the contribution from the appropriate root. One 
can show that the discarded root actually corresponds to 
the advanced Green's function, which implies that the 
boundary conditions, which were completely fixed by the 
choice of Legendre functions in the mode-sum expression, 
became ambiguous again in the process of summation. 

If we denote the appropriate root as V&O) we obtain the 
following closed-form expression for the Green's function 
for a static source in the Kerr black hole space-time: 

G sta tic(x,x') 

b 1 
~ 4^ A 1 AAV2 s inh(^|^o ~ A0|) - ab sin 9 sin 9' sin* 

(3.26) 

where we have reverted back to Boyer-Lindquist coordi- 
nates. 
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IV. THE SELF-FORCE ON A STATIC SCALAR 
CHARGE IN KERR SPACE-TIME 

In this section, we shall use our expression for the static 
Kerr Green's function obtained in the previous section to 
calculate the self-force on a static (with respect to an un- 
dragged, static observer at infinity) scalar particle in the 
Kerr background space-time. Since we have a closed- 
form expression for the appropriate Green's function, we 
can derive completely analytic expressions for the self- 
force. We restrict our attention to the simplest case of a 
point-like scalar charge q of mass m coupled to a mass- 
less scalar field $>(x). This calculation has previously 
been considered in [1] using the mode-sum regularization 
prescription (MSRP), developed by Barack, Ori and col- 
laborators [2"2"H2"5] , as well as Detweiler and collaborators 
The conjectured result for the self- force in [T] is 
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This derivation, however, is somewhat unsatisfactory in 
that the authors are forced to conjecture the unknown 
MSRP coefficient D a based on its known form in other 
space-times as well as inferring from numerical results 
that the contribution from the tail-integral is zero. By 
calculating the self-force by other more direct means, we 
can establish whether or not the authors' conjecture for 
the coefficient D a is correct (although this coefficient as 
well as higher order MSRP coefficients have been recently 
computed for the Schwarzschild [16 and Kerr [T7] space- 
times) and verify the accuracy of the numerical results. 
Moreover, since the MSRP appears to be an efficient 
method for calculating the self-force for more difficult 
trajectories such as circular and eccentric orbits in Kerr 
[5], it is important to establish the accuracy of this 
method by comparing with the results obtained by alter- 
native methods. More recently, a method known as the 
m-mode regularization prescription has been developed 
and applied to the Kerr case [291431] and again the calcu- 
lation presented in this section may serve as a standard 
non-trivial check of the accuracy of the m-mode regular- 
ization scheme. 



General Considerations 



The massless scalar field, $(x), satisfies 



(□ -£R)$(x) = -4tt P (x) 



(4.2) 



where □ is the wave-operator (d'Alembertian) on the 
background geometry, R is the Ricci scalar which is zero 
for vacuum space-times such as Kerr and £ is the cou- 
pling of the field to the curvature. The charge density, 
p(x), of the point-particle is 
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where z(t) describes the worldline 7 of the particle 
with proper time r, g = |det(<? h)| where g a b are the 
background metric coefficients and S 4 (-) is the four- 
dimensional Dirac distribution. The scalar field, &(x), 
moves on null geodesies of the background space-time 
whereas the scalar particle itself is massive and there- 
fore moves along time-like curves of the background. In 
addition, the field may scatter off the curvature of the 
space-time and re-interact with the particle. This radia- 
tion reaction process gives rise to a self-force 
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This self-force then appears on the right-hand side of the 
equations of motions for the scalar particle 

ma b = (g bc + u b u c )f s c cli = q(g bc + uV) V C $ R (z(r)) 

(4.5) 

where u b = dz b /dr and a b = Du b /dr are the four- 
velocity and four-acceleration of the particle, respec- 
tively. The crucial step in obtaining the correct self-force 
is identifying the correct radiative field that is regular 
at the particle's position, which we have called $r . The 
mass, m, appearing in the equations of motion is the 'dy- 
namical' particle mass, which in the scalar case evolves 
according to 
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An expression for the derivative of the radiative field 
that is regular at the particle's position can be obtained 
in terms of an integral of the retarded Green's function 
over the entire history of the particle's motion: 
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where R a b is the Ricci tensor of the background metric 
and a = Da b /dr is the Fermi derivative with respect to 
proper time of the four-acceleration. The last two terms 
are purely geometrical and easily evaluated. The global 
radiative term, the so-called tail integral term, $* ai1 , is 



$* ail (z(r)) = q lim 



W a G rct {z{T), Z {T'))dT' (4.8) 



where G TC t(z, z') is the retarded scalar Green's function 
satisfying 

(□ - ZR)G iet (z, z') = -g- l 'H\z - z'). (4.9) 

Since the Green's function is singular at coincidence 
z — z' , we require some limiting process in order to reg- 
ularize the self-force. When z is within a normal neigh- 
bourhood of z' , then the retarded Green's function has 
the following Hadamard representation [32] : 
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the difference between the retarded and singular fields: 
$ R (x) = $ ret (a;) - $ s (x) 

= 4nq / G rc t(a;, z(r))dr - 4irq / G s (x,z(t))(1t, 

J 'y J *y 

(4.11) 

where in order to calculate the self-force the radiative 
field must be evaluated on the world-line of the scalar 
particle. Detweiler and Whiting [33J have identified the 
singular Green's function that yields the correct radiative 
field which has the following Hadamard construction: 



FIG. 1: The retarded field at x depends on the entire par- 
ticle's history from — oo < r < T re t- The advanced field at x 
depends on the entire future of the particle from r a d v < t < oo. 
The singular field depends only on the trajectory of the particle 
between the advanced and retarded points. 



where a(z, z') is the Synge world function which is 
half the square of the geodesic distance between z and 
z', 5 + (<t) is the past light-cone delta distribution and 
0+(— c) is the past light-cone step function. The biscalar 
A(z, z') is the Van-Vleck Morrette determinant which is 
regular and symmetric, as is V(z,z'). The presence of 
the 5+(a(z, z')) in the direct part of the retarded Green's 
function implies that the direct part has support only on 
the past light-cone and therefore does not contribute to 
the self- force since the tail integral in Eq.(4.8l is entirely 



inside the past light-cone. Hence, it is the tail part of the 
Green's function that is responsible for the self-force. 

The causal regions are defined by the intersection of 
the past and future null cone of an arbitrary field point, 
x say, with the world-line, 7, of the scalar charge. Follow- 
ing Poisson [32] , we shall call the intersection of the past 
light cone of x with the world-line 7 the retarded point as- 
sociated with x and we denote it x' = z(r lct ), where T rot 
is the retarded time. Similarly, the future light cone inter- 
sects the world line at the advanced point x" = z(r a dv), 
where r a d v is the advanced time. The retarded Green's 
function depends on the entire history of the particle 
up to the retarded point. One can define an advanced 
Green's function in an analogous way where the depen- 
dence is now on the entire future of the particle from the 
advanced point to infinity. This dependence is seen in 
Fig0 



B. The Radiative Field for a Static Particle 



Rather than adopting the expression (4.7) for the gra- 
dient of the radiative field, in the static case it is more 
direct to begin with the definition of the radiative field as 



G s (x,z) = ^-\A 1/2 (x,z)6(a(x,z)) + V(x,z)Q(a(x,z)) 

The singular (direct) part of this Green's function has 
support only on the light cone due to the 8(a) term. For 
z(t) < z(r rct ) and z(t) > z(t^ v ), it is clear that x and 
z(t) are time-like related, meaning that a(x, z(t)) < in 
these regions. Moreover, for r rct < r < r ac i v , the points x 
and z(t) are space- like related and therefore a(x, z(t)) > 
in this region. These considerations combined with the 
presence of the Q(a) in the tail term of Eq.(4.12| imply 
that the singular Green's function has support only in 
the region r ret < r < r adv (see Fig[l]). 

We consider now the retarded field appearing in 



Eq.(4.11), we take the coordinates of the scalar charge 
at some proper-time r to be x a = z a (r) and change the 
integration variable to run over the coordinate time i us- 
ing u dr = dt, which gives 



$ ret (a;) = 47rg 



G ret (A*;x,x)£, (4.13) 



where At = t — t, x is the spatial position of the field 
point and x is the spatial position of the scalar charge. 
Generally speaking, the spatial points will depend im- 
plicitly on the time, x = x(i), but in the static case, by 
definition, there is no temporal dependence. For a sta- 
tionary metric, the Green's function may be decomposed 
into Fourier frequency modes: 



1 />oo pOO 

t (x) = Airq-^ / / e- !wAi G( 3 )(x, x) duj dt, 

J —OO J —OO 

(4.14) 



where Glj is a three-dimensional Green's function on 
a dimcnsionally reduced metric. For a static particle, 

(3) ~ 

Gu is independent of t and this integration may be per- 
formed, yielding a S(uj) term. The uj integral is now triv- 
ial yielding 



1 



$^f lc (x) = 4Trq-G { 3) (x,x) = 4^G static (x,x) 

(4.15) 

where, in the particular case of the Kerr geometry, 



1 



G sta tic(x,x) is given by Eq.(3.26) 
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Consider now the singular field. From Eq.(4.11) and 
Eq.(|4.12|), we obtain 



i r°o . 

<S>s(x) = -qJ {A 1 / 2 (x,x)5(a(x,x)) 

- V(x,x)Q(a(x,x))\dT. (4.16) 

In the first integral here, we change the variable of inte- 
gration to a using a- a u a dT = da, where by <7 ;a we mean 
covariant derivative with respect to the coordinates x a . 
The only contributions are when a = 0, i.e. at the ad- 
vanced and retarded points. The second integral van- 
ishes outside of the range enclosed by these two points as 
a consequence of the step function. Therefore, we have 



<f> s (x) 



^A 1 / 2 ^') t A^ 2 (x,x") 



2r rc 



2 ?"adv 



V(x, x(r))dr 
(4.17) 



where, for a static particle, r rct = cr.fU* is the re- 
tarded distance between x and the world-line and r a d v = 
—a t tnu is the advanced distance between x and the 
world-line. The integral term is regular in the limit as x 
approaches the world-line. The geometrical nature of the 
direct part is now explicit as the world-line is approached 
(r 0). 

Therefore, for a static particle in a stationary space- 
time, the self-force is given by 



/.s< 
a 



self/ 



q lim 



/4-7T 



Va(^Gi 3 2 ( X ,x) 



A 1 / 2 (x,x') 



2 TV 



A 1 ' 2 {x 1 x") 



2r 



adv 



V(x, x(r))dr 



(4.18) 



(3) 

where G^Iq is the zero frequency Fourier mode of the 
four-dimensional retarded Green's function. 



C. Calculating the Retarded and Advanced 
Distance for a Static Particle 



In this Section, we shall outline our method for cal- 
culating the retarded and advanced distances, r rot and 
''adv, respectively, for a static particle in a stationary 
space-time. Calculating these quantities requires eval- 
uations of a derivative of the geodesic distance at the 
retarded point x' and the advanced point x" , which are 
connected to the field point a; by a null geodesic defined 
by a(x,x) = 0. One may obtain a coordinate expansion 
of a by assuming the form 

2a(x, x) =g ab Ax a Ax b + A abc Ax a Ax b Ax c 

+ B abcd Ax a Ax b Ax c Ax c + 0(Ax 5 ) (4.19) 



and then substituting this ansatz into the defining equa- 
tion for a 



2o = g ab o, a o, b . 



(4.20) 



Equating powers of Ax yields 



A nb r — 



1 



* 9(ab,c) , 



B n hrd — 



1 



abed ^9{ab,cd) 



l 
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nun 



9efL (ab 



(4.21) 



where, in our conventions Ax a — x a — x a , i.e., field point 
expanded about the particle world-line, and round brack- 
ets around indices imply symmetrization. Then 

2cr t - =2g tt At + 2g at Ax a + 3A ttt At 2 + 6A att AtAx a 

+ 3A apt Ax a Ax f) + 4B tttt At 3 + \2B attt Ax a At 2 
+ l2B am Ax a Ax l3 At + <iB aMt Ax a Ax f3 Ax 1 
+ 0{Ax 4 ) (4.22) 

which we wish to evaluate at the retarded and advanced 
points. This requires two additional steps: First, the 
tensors g ab , A abc and B abc( i are each evaluated at the 
field point x. However, we wish to evaluate them at the 
retarded and advanced points which are on the parti- 
cle's world-line. Therefore, we Taylor expand these coef- 
ficients about the world-line. Second, we solve the equa- 
tion a{x,x) = iteratively to obtain the advanced and 
retarded times as expansions in Ax a . Substituting into 
our expressions for a j gives us the advanced and retarded 
distances. 



Taylor expanding the coefficients in expression (4.22) 
about the world-line yields 

2cr t - = e(2g tt At + 2g at Ax a ) + e 2 {-2g ttta Ax a At 

- 2g at ^Ax a Ax 13 + SAutAt 2 + 6A att AtAx a 
+ 3A a ptAx a Ax^) + e 3 (g tt , al3 AtAx a AxP 

+ g at ^Ax a Ax fi Ax~< - ?>A tU a Ax a At 2 

- QA att . B AtAx a Ax p - ZA afitn Ax a Ax p Ax"* 



+ 4B tttt At 3 + \2B aUt At 2 Ax a 

+ l2B a0tt Ax a Ax fi At + 4B aMt Ax a Ax? As 7 ) + 0(e 4 ), 

(4.23) 

where an overbar means that the tensor is evaluated at 
x which is on the world-line and we have introduced a 
book-keeping parameter, e, in order to keep account of 
powers of Ax. It is worth noting that there is no need 
to distinguish metric dependent tensors g ab , A abc and 
B a bcd from those evaluated at x' or x" since these tensors 
depend only on the spatial coordinates for any metrics of 
interest, which have the same spatial dependence for a 
static particle (see Fig [2| . The difference in evaluating 
Eq.(4.23) at various points on the world-line of a static 



particle arises in the At term. 

The difference At corresponding to the advanced and 
retarded times we shall denote Ai a( j v and At rct , respec- 
tively. These times are defined by the intersection of the 
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Upon substitution of Eq.(4.25) into (4.24), we obtain 



(t,x) 



(*,*) 



FIG. 2: The spatial coordinates of the world-line of a static 
particle are fixed over coordinate time. 



2a(x,x) = e 2 {g tt At\ + 2g ta At 1 Ax a + g a0 Ax a Ax p ) 
+ e 3 {A ttt At\ + [3A tta ~ gtt, a }AtlAx a 
+ [3A ta p - 2g ta)j9 ]At 1 Ax a Ax^ 
+ [A a01 - g aP ^]Ax a Ax f3 Ax' 1 + 2g ft At 1 At 2 
+ 2g ta At 2 Ax a ) + e 4 (B tttt Ati 
+ [4B ttta - A m JAt\Ax a 
+ fiBttaf} ~ SAtta,? + \gtt. al 3\AtlAx a Ax 13 
+ [4B taPl - 3A ta/3 ^ + gt a ,p y ]Ati Ax a Ax fs Ax 1 
+ [B a0lS - A a M,& + 2 -g a p n s]Ax a Ax p Ax~ 1 Ax s 
+ 3A ttt At 2 At 2 + [QA tta - 2g tt . a ]At 1 At 2 Ax a 
+ [3A taP - 2g ta . p ]At 2 Ax a Ax 13 + 2^^^ 
+ g tt At\ + 2g ta At 3 Ax a ) + 0(e 5 ). (4.26) 



null light cones with the world-line 7 and we are there- 
fore required to solve a(x,x) = as a function of At. A 



coordinate expansion for a has been given in Eq.(4.19) 
along with expressions (4.211 
pendent tensors A a b c 



As before, the metric dc- 
and B a b c d ought to be evaluated 



on the world-line and so a Taylor expansion about the 
world-line is required, yielding 



Setting a(x,x) equal to zero and equating equal orders 
of e gives an iterative solution for At. In particular, at 
0(e 2 ), we have 



g tt At( + 2g ta At 1 Ax a + g af3 Ax a Ax fs = 0, (4.27) 



2a(x, x) = e 2 (g tt At 2 + 2g ta AtAx a + g af3 Ax a Ax 13 ) 
+ e 3 (A ttt At 3 + [3A tta - g ttta }At 2 Ax a 



[3Atat3 

[A a p~f - 



2g ta ,p}AtAx a A 



X' 



g al3 , 1 }Ax a Ax p Ax 1 ) 
f4 1 A tttia ]At 3 Ax a 



[6B«a/J - 3A tta ,f3 + \g ttMi \At 2 Ax a Ax p 
[AB taM - 3A taf3n + g ta ^]At Ax a Ax@ Ax 1 

[B af)l s - A a p-y,8 + \g af }. 1 s]Ax a Ax p Ax 1 Ax s ) + 0(e 5 

(4.24) 



which is a quadratic in Ati with solutions 



At* 



± _ -g ta Ax a ± y/fapAxP) 2 



9tt9oi0 



Ax a Axf~ z 



9tt 



(4.28) 

These two solutions correspond to the first order approx- 
imation to the advanced and retarded times. From Fig [2] 
it is clear that At rot = (t' — t) <0 and At adv = (t" - 1) > 
0. Since <jta < (a = 0...3), the appropriate identifica- 
tion of the retarded and advanced times is 



Atl.ret — At^~ 

We wish to solve a = iteratively for At so we further At% a d v = At^. (4.29) 

expand 



At = Ati + eAt 2 + e 2 At 3 . 



The higher order corrections are linear in the dependent 
(4.25) variables and so it is straight forward to calculate the 
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0(e 3 ) 
At 2 . 



and 0(e 4 ) solutions: 
1 



At 



2(s„Ati, ret + g ta Ax a ) 
[SA tta - g tt , a }Atl Ict Ax a 
[3A taP - 2g ta ^]At 1}Iet Ax a Ax p 
[A aM - g af)n }Ax a Ax p Ax~< 
1 



(4.30) 



2.adv 



2(ff tt At 1Mv + g ta Ax 
[3A tta - g t t,a]Atl adv Ax a 
[3A ta p - 2g ta ^]At 1Mv Ax a Ax 

[A aM - g al 3, 7 }Ax a Ax p Ax~ ( 



^MtuAtl 



adv 



At, 



(4.31) 



2[g tt At ltTet + g ta Ax a ) V 
[4Bttta " JA^Ax" 
[6B Mq/3 - 3A tta ,p + \g t t, a p]Atl , ct Ax a Ax fi 
[4£? to| g 7 - 3A tQ(3!7 + g tQ ,/3 7 ]Aii !rot Aa^Aa^Ax 7 
[B aM5 - A aM j + y al3 ^ s ]Ax a Ax^Ax^Ax s 
[6A Ua - 2g t t, a }At lA . ct At 2 ,rctAx a 
[3A tafj - 2g ta ^]At 2 ^ t Ax a Ax 13 



3A ttt At 2 ret At 2) ret + 9tt A*Lc 



(4.32) 



Ai 



1 



3, adv 



adv 



2(g tt At 1Mv + g ta Ax° 
[4B«ta - ^.JA^ adv Ax a 
[6fl«a/3 - 3vW/3 + ig tt , a/ 3]At2 adv Aa; Q Aa ; /3 
[4B ta/ 3 7 - 3A ta n + g ta ,fs-y}At 1Mv Ax a Ax^ Ax 1 
[B aMS ~ A aPl ,s + lg a p, 7 s]^x a Ax^Ax^Ax 5 
[6A tta - 2g tt , a ]At 1Mv At 2Mv Ax a 
[3A ta p - 2g ta ^]At 2Mv Ax a Ax f) 



3A t «At? iJldv Ata 1 ret + 9tt At| 



adv 



(4.33) 



We are now in a position to write down expressions 
for the retarded and advanced distances as coordinate 



expansions. Combining Eq.(4.23l and Eq.(4.25) with 
Eqs. (14.281)- (14.331) implies 



1 



^ret = CTt'U 



2yf=mt 



(4.34) 



■Rri 



AB tttt Atl Ict + [\2B ttta 



3A tt t, a }Ati tTet Ax a 
+ [UBttop - 6A tta ,p +gtt, a p}At 1 . rct Ax a Ax> 3 
+ [4B taPl - 3A taPn + g ta ^]Ax a Ax' 3 Ax' y 
+ 6A ttt At hl . ct At 2:Ict + [6A tta - 2g tt:a ]At 2trct Ax a 
+ 2g tt At 3 , ret . (4.37) 

Analogously 



^*adv 



.t" 



2V = 9u 



-^adv^ ~t~ Qadv^ 



^adv£ 3 

(4.38) 



where P a dv, Qadv and i? a dv are obtained from Eqs. (4. 35) 



(4.37) by replacing the Ai rGt expressions with At adv 



terms. 

We may now use these expressions for the advanced 
and retarded distances in Eq.(4.18). We first note that 
covariant expansions for the bitensors A 1 ^ 2 (x,x) and 
V(x,x) are found in |34j . which for a massless field in 
a Ricci-flat space-time are given by 



A 1 ' 2 {x 1 x) = \ + 0{Ax A ) 
V[x,x) = 0{Ax 2 ), 



(4.39) 



which implies that the integral term involving the biten- 
sor V(x, x) in Eq.(4.18) will vanish in the limit x — > x 



and so will not contribute to the self-force. The self-force 
reduces to 



/r lf = q 2 lim 



4n G (3) fxx) 1 1 

U Z 7*rct * ^*adv , 

(4.40) 

Substituting our expansions for the retarded and ad- 
vanced distances in Eq.(4.34| and Eq.(4.38), respectively, 



we obtain the self-force on a static, massless scalar charge 
in a stationary, Ricci-flat space-time: 



ff { = W lim 



V-gtt V a ( G^ 3 2 (x, x) - G sing (x, x) 



(4.41) 



where 



l r l 



47T 

(Q 



l 



Qrct . Q 



adv 



-^ret fadv 
t — ^ret-Sret) 



P 2 



pi 
^adv 



(Qadv -Padv-Radv) 



p3 
■'ret 



p3 
adv 



(4.42) 



It is easy to check that Eq. ( 4.41 )-( 4.42 ) gives us zero 



self-force for a static particle in Schwarzschild, as ex- 
pected ng . 



where 



P rct = 2g tt At hrct + 2g ta Ax a , (4.35) 
Qret = 3A ut Atj rct + [6A tta - 2g tt , a ]At hl . ct Ax a 
+ [3A taP - 2gta.js} Ax a Ax^ 1 + 2g tt At 2 ^ ct , (4.36) 



D. Self-Force on a Static Particle in the Kerr 
Space-Time 

We now calculate the self-force on a scalar charge with 
arbitrary fixed spatial Boyer Lindquist coordinates out- 
side the ergosphere (the notion of a static particle breaks 
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down on the ergosphere where the Killing vector d/dt 
becomes null). The algebra involved is unnecessarily 
cumbersome for arbitrary separations, so we calculate 
each component of the self-force separately, which allows 
us to set two of the three separations to zero. For ex- 
ample, if we wish to calculate the radial component of 
the self-force, then from Eq.(4.41) we require partial ra- 



dial derivatives of the Green's function and the singu- 
lar field before taking a coincidence limit. However, in 
this case, one can take partial coincidence limits Acp = 
and A9 — from the outset since these terms just pass 
through radial derivatives. This procedure makes the cal- 
culation more palatable from an algebraic perspective. 

There are several steps involved in the calculation: 
First, we derive analytic expressions for \l/o for small 
separations. We then use our results for to obtain 
analytic expressions for the Green's function for the var- 
ious separations. Then the singular terms are calculated. 
Finally we compute derivatives of the regular field and 
take coincidence limits to obtain the self-force. 



1. Calculating $o for Small Separations 

In order to calculate the Green's function up to order 
O(Ax) (which is what is required for the self- force), we 
require "Po up to 0(Air 3 ), where V^o satisfies the tran- 
scendental equation ( |3.18| , which in Boyer Lindquist co- 
ordinates is 



cosh f-Oo - A<t>) 
\a 

where 



(r — M){f — M) — b 2 cos A 



A l/2 A l/2 



cos A = cos 9 cos 9 + sin 9 sin 9 cos ■ 



(4.43) 



(4.44) 



As discussed in the previous section, we cannot solve this 
equation analytically in general. However, we shall show 
that it can be solved for small separations up to any 
order. Reintroducing the book-keeping parameter e to 
keep track of orders of Ax, we write 



4# 3 e + ^ 2) e 2 + ^ 3) e 3 + 0(e 4 ). (4.45) 



Substituting this into Eq.(4.43), writing r — r — eAr, 



9 = 9 — eA9, <j> = (f> — eAcft, expanding about e = 
and then equating equal orders of e gives a simple set 
of recursion relations for Vfg , ^o 2 ' e ^ c - For example, 
equating powers of e 2 yields the following quadratic in 



(1)>|2 



where 



2A* 



(i) 



A(f> 



A(A - a 2 sin 2 9) K ' 



Ar l + A9 2 A. 



(4.47) 



The two solutions of this quadratic are 



(i) _ 
o,± — 



A Aqb j a 2 Asm 2 9 AcfP a 2 s 2 

A -a 2 sin 2 9 V (A - a 2 sin 2 8) 2 + A(A - a 2 sin 2 0) ' 

(4.48) 

where we recall that is the particular choice that 
yields the correct singular field, and as explained in 
Sec. [Hi] corresponds to choosing the more positive so- 
lution: 



(i) 



A Aqb 



a 2 Asin 2 0A0 2 



A - a 2 sin 2 9 ]j (A - a 2 sin 2 8) 2 A(A - a 2 sin 2 9) ' 

(4.49) 

The higher order coefficients are linear in the depen- 
dent variable and therefore uniquely defined in terms of 
the appropriate solution. However, for arbitrary 

separations they become increasingly cumbersome from 
an algebraic perspective so we calculate each component 
of the self-force individually by specializing to a particu- 
lar separation. 

Throughout the remainder of this section we shall as- 
sume, without loss of generality, that Ar > 0, A9 > 0, 
Acp > and a > 0. Then the limit in Eq.d4.41h is a 



one-sided limit from the left. In fact, the two-sided limit 
does not exist since the Green's function expansion for 
small separations involves terms like 
For radial separation, we obtain 



*«(Ar) 

^ 2) (Ar) 
^4 3) (Ar) 



aAr 



AV2(A-a 2 sin 2 9) 1 / 2 

a(f - M)(2A - a 2 sin 2 9)Aj 



2A3/2(A-a 2 sin 2 (9) 3 / 2 



"24A 5 / 2 (A-a 2 sin 2 9f' 2 
-a 2 A sin 2 9(m 2 + a 2 + 



3a 4 sin 4 6>(4b 2 + 3A) 



4A 2 (6A + 86 2 ) 
(4.50) 



-24A) 



which along with Eq. ( 4.45 ) defines for radial separa- 
tion up to 0(e 3 ) terms. 

Similarly, for 0-separation, we have 

>(!)/•- aM 



n 1 m = 

¥*\A9) = 
^ q 3) (A9) = - 



(A-a 2 sin 2 0)1/2 

a 3 cos 9 sin 9A9 2 
~ 2(A - a 2 sin 2 Of/ 2 



A{A + b 2 ) 



24(A - a 2 sin 2 9) 5 / 2 
4a 2 sin 2 (9 ( A ~ 2a 2 )+ 4a 4 sin 4 # (4.51) 
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while for ^-separation, we have 



A0A 1 / 2 



(A 1 / 2 _ a s i n | 



^ 2) (A0) = O 

(3) aA^sin^A + frW^ 3 
^ (A ^ = Mp^kf ' (452) 



small separations up to O(Ax) (or 0(e)). 

We begin with radial separation. The Green's function 
expression Eq.(3.26) for radial separation is given by 



Gl mic (r, r, 9) 



1 



4tt AV2AV2 sinh[f * ] - ab sin 2 9 sin * c 



1.53) 



2. The Green's Function for Small Separations 

We now employ these expansions for m order to 
calculate the retarded Green's function expansions for 



where * = (Ar)e + (Ar)e 2 rf' (Ar)e 3 + 0(e 4 ) 



,(2) 



(3), 



with <3> 



(i) 



etc. given by Eq.(4.50). Writing r = 



r — eAr and doing a series expansion about e = (this is 
equivalent to a small Ax expansion), we obtain 



Cstatic( r > f, 



47T 



+ 



A 1 / 2 



a 2 (f - M) sin 2 



(A - a 2 sin 2 9y/ 2 Ar 2A 1 / 2 (A - a 2 sin 2 9f/ 2 
2 sin 2 9[A{AA + Ab 2 - a 2 ) - a 2 sin 2 9{A + 2b 2 )} Ar 



8 A3/2(A- a 2 sin 2 9f> 2 



0(Ar 2 ) 



where we have set e = 1. 



For ^-separation, Eq.(3.26) becomes 



^static( r i ^' 



47r Asinh(|* ) - absm{9 - eA9) sin sin \& 



where is given by Eq.( 4.45 (and (4.51). Expanding about e = gives 

1 



G" static (f,9,9) 



1 

47T 



a' cos V sin t 



(A-a 2 sin 2 0) 1 /2A0 2(A - a 2 sin 2 Of/ 2 
[A (A - 3b 2 ) + 2a 2 sin 2 0(3a 2 - 4A) - 2a 4 sin 2 0] A0 



24(A - a 2 shr 



2M5/2 



O(A0 2 ) 



For ^-separation, the Green's function is 

Gl tic M,A0) 



1 



4tt Asinh[^(* o -A0)] — a 6 sin 2 sin ^ 
where is given by Eq.(4.45) and (4.52). In this case, the series expansion for the Green's function is 

1 A<j) 



G*(r,9,A4>) 



1 

47T 



A 1 /2 s i n 0A0 24A 1 /2(A-a 2 sin 2 6») 3 
3A 2 sin 9(M 2 + a 2 ) - 8a M 2 A 3/2 sin 2 9 



A* csct 



3a A sin" 



0(2M 2 - a 2 ) + a 4 b 2 sin 5 #} + 0(A</> 2 ) 



(4.54) 



(4.55) 



(4.56) 



(4.57) 



(4.58) 



3. Analytic Expression for the Self-Force 

We now have expansions for both the Green's function and the singular field required to compute the self-force. 
As mentioned above, a particular component of the self-force may be evaluated by taking partial coincidence limits 
inside the derivative. For example, the radial component is given by 



iV sclf = 4tto 2 lim 



V^tht d r (G r static (r, r, 9) - G r sing (r, r, 9)) 



(4.59) 
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where Gg tatic (r, f , 0) and Gg ing (r, f, 0) are the radially separated Green's function and singular field, respectively. 
We have computed the Green's function for the various separations in the previous section, the singular field for a 
particular separation is given by Eq.(4.42) in the partial coincidence limit. After some algebra, we obtain 



G r sin Jr,f,6) 



r e (r I 



Gf ing (f,0~,A0) 



1 

47T 
+ 

1 

47T 

+ 

1 

47T 



A 1 / 2 a 2 (f-M) sin 2 

(A - a 2 sin 2 0)V2 Ar + 2AV2(A-a 2 sin 2 0) 3 / 2 
a 2 sin 2 0[A(4A + 46 2 - a 2 ) - a 2 sin 2 0(A + 2b 2 )] Ar 
8A 3 / 2 (A-a 2 sin 2 0) 5 / 2 
1 a 2 cos 9 sin 

(A - a 2 sin 2 0)V 2 A0 ~ 2(A - a 2 sin 2 0) 3 / 2 
A(A - 36 2 ) + 2a 2 sin 2 0(3a 2 - 4A) - 2a 4 sin 2 0] A0 
24(A - a 2 sin 2 Of/ 2 
1 Ac/) 
A 1 / 2 sine Ac/) + 24A!/ 2 (A-a 2 sin 



2 ms 



3A 2 sin 0(M 



2 1 a 2 ) 



3a 2 Asin 3 (9(2M 2 



A 3 esc 



a 4 6 2 sin 



0(Ar 2 ) 



O(A0 2 



A + O(A0 2 



(4.60) 



(4.61) 



(4.62) 



Gf 



3 (r,f,0")-GI inK (r,f,0") = O(Ar 2 ) 



and therefore Eq.(4.41) clearly gives 
/f lf = 0. 



Comparison of Eq.( |4.54 1 with Eq.(4.60), we see that 

(4.63) 

(4.64) 

(4.65) 

(4.66) 



Similarly, from Eqs.(4.56l and ( 4.61| ), we have 
G s 9 inK (f,0,0) = O(A0 2 ) 



and hence Eq.(4.41) yields 

/r = o. 



Therefore the only component of the self-force lies in the 
azimuthal direction, as suggested by the numerical re- 
sults of Burko and Liu [T]. Subtracting Eq.(4.62) from 
Eq.(4.58) gives the regular field 



G 



a , 



Gl 



static V ' " 5 TV " sing 

M 2 A sin 2 



12tt (A 



a 2 sin 2 0) 3 



(f,0,A^) 
0(A^ 2 ) 



(4.67) 



Differentiating with respect to <fr, including the 4 71-/^ /— g tt 
term and taking the limit <f> — > <j> according to Eq.( 4.41| 



gives the azimuthal component of the self-force on a static 
particle in the Kerr space-time: 



M 2 



a 2 sin 2 0W2£i/2 



(4.68) 



Our results clearly agree with the numerical results of [T] . 
We infer from this agreement that the MSRP parameter 
D a that the authors of Ref. [I] conjecture is correct, at 
least for a static particle. Our analytic results ought 
to offer a good standard test for any numerical regu- 
larization prescription in the non-spherically symmetric 



case. We note also that the result presented here is triv- 
ially extended to the Kerr-Newman black hole by taking 
M 2 -> M 2 - Q 2 , where Q is the charge of the black hole. 



Combining the results Eqs.(4.64|, (4.66) and (4.68) 



and assuming, for notational convenience, that the scalar 
charge is located at Boyer-Lindquist coordinates (r, 9, <f>) 
gives the following representation of the self-force on a 
static scalar charge in the Kerr background geometry: 



/ b sclf = ^ 2 aAsin 2 , 



, HP 

(A - a 2 sh?0j572si72 



8\. (4.69) 



CONCLUSIONS 



In this paper, we derived a completely closed-form rep- 
resentation of the retarded Green's function for a scalar 
particle at fixed spatial Boyer-Lindquist coordinates in 
the Kerr space-time. This representation relied on two re- 
sults for the Associated Legendre functions, one of which 
is derived in the Appendix and the other we derive else- 
where ) 1 3 j - These formulae will likely prove useful in 
other contexts. As a check of the validity of our closed- 
form Green's function, we computed, analytically, the 
self-force on a static scalar particle in Kerr space-time 
and showed that our result agrees with the conjectured 
result of Burko and Liu [1] . The calculation presented in 
this paper is therefore the first proof of the result implied 
by the numerical work of Ref. [T] . 
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Appendix 

Summation Formula for the Product of Associated 
Legendre Functions 

In this appendix, we shall derive the following summa- 
tion formula for Associated Legendre functions: 

HH {l + m) - 

2 cos(mcos _1 (C)) , , 

~ ^(sin^sin 2 6" - (x - cos 9 cos 9') 2 ) 1 / 2 ^ ' 

where 

x — cos 9 cos 9' 
= sin 9 sin 0' ^ ' 

and x — cos A must lie in the range 9 — 9' < A < 
min{0 + 9',2tt — 9 — 9'}. For A outside of this range, 



the series sums to zero. Though not a new result, it is 
not well known, the only reference of which we are aware 
that includes this formula is Ref. [TH]- We have there- 
fore included a proof of this result, following closely the 
method of Baranov [2D] who has proven the m = case, 
as the technique might be valuably extended in other 
directions. We begin with the well-known Christoffel- 
Darboux formula [36] 

_ n + 1 P n+1 (x)P n (t) - P n (x)P n+1 {t) 

~ 2 x-t ■ W 

Taking 

x = cos A, 

t — cos [i = cos 9 cos 9' + sin 9 sin 9' cos <p, (4) 
we may employ the Legendre addition theorem to obtain 



2 ( / I 7T& } . 2 /■ 

;=0 m=-l v ' 

Reversing the order of summation, multiplying across by cos(m'<^) and integrating from to 2tt gives 



7T V(2/ + l) 7r - ^yP; (cos6»)P ; (cos(9 )P ; (x) = — - — / cos(m^) dy>. (6) 

L — \ra\ 

We change the integration parameter from tp to t using 

*P = ± -> dt = (?) 

v/(cos(6> - W) - t)(t - cos(6» + 9')) 

where the choice of sign here is determined by the fact that t, as a function of <p, is symmetrical about (p = tt, 
decreasing from cos(9 — 9') at cp = to cos(# + 9') at <p = w, and then increasing symmetrically to cos(9 — 9') at 
ip = 2n. Hence, we obtain 

77 E (2l + l)f—^Pr(cos9)Pr(cos9')P l (x) 

l=\m\ ^ >' 

,c«(fl-e') P n+1 (.T)P„(t)-P»(:r)P» +1 ft) 

= m+l / cos togs it)) dt 8) 

Ja»(p+o') (x - t)^(cos(9 - 9') - t)(t - cos(0 + 9')) 



where 



, ft — cos 9 cos 9'\ 

cp(t =cos- x , , . 9 

V sm sm 6r / 

In order to obtain the result (JlJ , we require an asymptotic analysis of the integral in Eq. ([8]) for large n so that the 
limit as n — > oo may be taken. Such an asymptotic analysis is complicated further by the fact that one needs to allow 
for the possibility of small or vanishing x — t (which only occurs when x lies inside the integration range). Using the 
large n asymptotic formula [14J 



1/2 TT / 2 \3/2 



P„(cosz)=( — ) C0S [(n+l>-f]-^( — ) cos[(n+|)z+f] + 0(n- 5 / 2 ), (10) 

Vmsinz/ lo\7msinz/ 



15 



we can deduce 



„ ,,„,,-, „ ( „ 2[sin(i(A-M))cos((n+l)(A + M))-sin(i(A + M))sin((n+l)(A-M))] , . 

P n+1 (x)P n {t) - P n (x)P n+1 (t) = — . f == V A n (x,t) 

ir(n + 1) ysm /x sin A 



(11) 



where t = cos fi, x = cos A and the remainder A n (x, t) is given by 

1 /cos[(n+|)A- f]cos[(n+ |)/x+f] cos[(n + |)A - f ] cos[(n + + f ] 



A n (x, t) 



47m 2 



y 7 sin A sin 3 fi 



\J sin A sin 3 fi 

cos[(n + |)A + f ] cos[(n + §)/j - f ] cos[(n + |)A + f ] cos[(n + - f ] 



0(n" 



(12) 



V sin 3 A sin fi y 7 sin 3 A sin fi 

One can choose a function of cc such that the trigonometric terms in the brackets are bounded for all values of t, i.e 



(13) 



where c(x) is independent of n and t. The approximation (111 with Eq.(13) is valid whenever x is not close to t but 
as mentioned above, we require an alternative estimate to that of Eq.(13) that accounts for small or vanishing x — t. 
To achieve this, we can use the asymptotic formula for the associated Legendre function of order 1 

if (cob z) = ( — ) V2 cos[(n + \)z + f ] + /_2n_\ 3 / 2 cos[(n + \)z + f ] + O^ 3 / 2 ), (14) 
Wsinz/ Ion 2 \7rsmz/ 



as well as Eq.(10l to obtain the estimate 

2[sin(|(A - fi)) sin((n + 1)(A + fi)) - sin(|(A + /*)) cos((n + 1)(A - fi))] 



P n (x)P^(t) - P n+1 (x)P^(t) 



n^/ sin fi sin A 



B n (x,t), 
(15) 



where 



B n (x,t) 



Given that 



1 /3cos[(n + |)A-f]cos[(n+|)/i+ 3 f] 3cos[(n + |)A - f ] cos[(n + §)/j + 



Ann 



3tt_] 
4 J 



y sin A sin 3 /i 



y 7 sin A sin 3 u 



cos[(n + |)A + f ] cos[(n + \)fi + f ] cos[(n + |)A + f ] cos[(n + + f ] 
y 7 sin 3 A sin /x y 7 sin 3 A 



sin /x 



0(n 



Pi it) 



sm u - 



d P(t) 
dt ' 



(16) 



(17) 



we may divide both sides of Eq.( 15 ) by sin fi and integrate from x to t to get 

2[sin(±(A - fi)) cos((n + 1)(A + fi)) - sin(|(A + fi)) sin((n + 1)(A - /i))] 



P„+l(aOP„(t)-P„(V)P n+ i(t) 



7r(n + l)ysin fi sin A 



* B n (x,t) 



dt. 



sm /x 



Comparing with ( 11 ), we see that 



A n (x, t) 



1 B n {x,t) 



dt. 



sm fi 



(18) 



(19) 



Again the trigonometric terms in B n (x,t)/ sin fi are bounded. If we take m to be the lower bound and M the upper 
bound and define the positive constant c\ — | max{m, M}\, then we get the following bound for A n 



\A n (x,t)\ 



sm fi 



< 



(20) 
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The geometric mean of Eq.(13) and (20) yields a new estimate for the remainder term: 

y/c{x)Cx\x - t\ 



\A n {x,t)\ < 



,3/2 



(21) 



Employing Eq.(|ll| with the estimate (21 1, we have 



cos(e-e') 



2 cos(m y (t)) 



t(e+fl') tt(^ — t)\/sin^sin A^/(cos(0 — 9') — t){t — cos(0 + 0')) 

sm(|(A - n)) cos((n + 1)(A + /x)) - sin(|(A + (i)) sin((n + 1)(A - //))] eft + 0(-=). 



It is convenient to integrate over /i rather than t, and express the singularity as 

11 ' /(A,/*) 



x — t cos A — cos fi (fi — A) sin A 



(22) 



(23) 



where / is a continuous function. Since / contains no singularities, and in view of Eq.(13), it is easy to see that the 
contribution of the integral containing / in the integrand is 0(n _1 ). Therefore, 



it {21 + 1 



{l-m)\ 



Pr{cos9)Pr(cos9')Pi(x) = 



l=\m\ 



(I + m)\ 

2 cos(m ip(fi)) 



sin /i 



d^ + 0{^=). 



(24) 



lfj,_ Try/ (cosp — W) — cos ^)(cos^i — cos(0 + 9')) V sin A 

rsin(|(A - fi)) cos((n + 1)(A + fi)) - sin(|(A + //)) sin((n + 1)(A - /*))- 

— A - 

where =9 — 9' and /i + = min{0 + 0', 2tt — 9 — 9'}. Without loss of generality, in what follows, we shall assume 
that 9 + 9' < 7T so that n + = 9 + 9' . 

The first integral in Eq. ( 24 1 may be written as 



cos(m 



sin fi d[i 



2 sin(|(A - (J,)) cos((n + 1)(A + fi)) 
ny/(cos(9 - W) - cosfi)(cosfi- cos(0 + 9')) V sin 3 A \i - A 



where q m is independent of n but not of m. One can further write 

A R 

9m (M) = 



q m (fj,) cos((n + 1)(A + 

(25) 
(26) 



V/i - (0 - 0') V^ + O'- A* 

for particular constants A m and B m and some continuous function q m> x (jj,). Considering the first singular term here: 
- cos((n + 1)(A + M )) ^ =2 f ^ cog((n + (A + e _ e , + g2))ds 



--2cos[(n + l)(\ + 9-9')} / cos[(n + l)s 2 ]ds 



/28' 



Vv e ') 



(27) 

where we adopted the transformation s = Wji — (9 — 9') in the first line. We further make the transformation 
S = s>/n + 1 which gives 



2sin[(n + l)(A + 0-0')] / sin[(n + l)s 2 ]ds 



/29' 



cos((n + l)(A + At )) _2cos[(n + l)(A + 0-0Q] fV^+^ 



Vn + 1 

2sin[(n+l)(A + 0-0')] / > ^ 2(n+1)e ' 
Vn + 1 



cos(S 2 )dS 
sm(S 2 )dS. 



(28) 
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These last two integral are the well-known Fresnel integrals and they converge for large values of the argument 
(n — > oo). Since cos(nir) / \fn — > and sm(nx)/y/n — > as n — ¥ oo, we must have 



cos((n + l)(A + /i)) 



^_ - (0 - V) 

An identical analysis for the other singular term in Eq.( |26[ ) reveals 

cos((n + l)(A + /i)) 



rf/i — > as n — > oo. 



V0 + 0'-v 



d/i =— > as n — > oo. 



(29) 



(30) 



Hence, the contributions from the singular terms in Eq.(26 



vanish for large n. 



Turning our attention now to t he r egular part of Eq.(26), taking our integral over the entire real line (since 



consideration of Eqs.(|8j), (11) and (13) show that q m ,i{n) is 0(n x ) outside the range of integration and therefore 
does not contribute in the large n limit), we obtain 

1 f°° If 00 
9m,i(A i )cos((n + l)(A + ^i))d/i = - / q m ,i((J,) cos((n + 1)(A + n))dfx + - I q m ^{v)cos{{n + l)(\ + v))dv. 

* J — oo ^ J — oo 

(31) 

" = !'-- ; , (32) 



Making the substitution 

in the second integral yields 

Qm,i (/■*) cos((n + 1)(A + fi))dfx 



n + 1 



QWiCaO - ImMH 77) cos((n + 1)(A + fj))dfJL (33) 

V n + 1 / J 



which vanishes as 



Finally, we consider the second integral in Eq.(24). This integral may be written as 
"+ sin((n +l)(/i -A)) 



fi — A 



-T m (fi)dfi 



M+ da 
{sin((n + l)(/i - A))T m (A) + sin((n + 1)0*- ~ T m (X)]} ' 



where 



cos(m <p(n)) sin(|(A + fi) 



sm/i 



7T v /(cos(6» - 9') - cos Ai)(cos Ji - cos(0 + 9')) V sin 3 A ' 
For the second integral in the expression above, one can write 

T m (fi) — T m (A) C m D m 



jJL — A 



y/fi-(0- 9') Vd + d 1 - A* 



■Pto(m) 



/x — A 
(34) 



(35) 



(36) 



where p m is a continuous function. Applying a similar analysis as that above, one can show that the contribution 
coming from this integral converges to zero as n — > oo. 



The only contribution to the sum in ( 24 ) therefore comes from the first integral on the right-hand side of ( 34 ) : 

,(„+l)(0+0'-A) • 



T m {\) 



sin((ra + l)(/i - A)) 



dfi = T m (A) 



/i — A 

In the limit as n — ¥ oo, this integral converges to the Dirichlet integral 

»(n+i)(9+e'-A) 



-dr. 



(37) 



(n+l)(6>-0'-A) 



lim 



T m (A) 



smr 



(n+l)(0-6'-A) 



= T m (X) 



SUIT 



dr = TiT m {\) 



2cos(m(p(X)) 



sin 2 9 sin 2 0' - (cos A - cos 9 cos 9') 2 



(38) 
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where 



ip(X) = cos 1 ( 



_ 1 / cos A — cos cos 9' 



sin v sin f 



(39) 



We note here that for A outside the range (/i_,/i+), the integral (34 1 would have no singularities and the integral 
would have vanished in the limit n — > oo. Therefore we have shown 



J=|m| 



(Z + m)! 



7T (sin 2 sin 2 0' - (x - cos 6> cos 6>') 2 ) 1/2 



(40) 



where 



c = 



X — COS tt COS f 



sin v sin I 



(41) 



and x = cos A must lie in the range 9 — 0' < A < min{6> + 0\ 2n — — 0'}. 
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